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Department of Physics, Case Western Reserve University, Cleveland, OH 44^06-7079 

We re-analyze constraints on the cosmological constant that can be obtained by 
examining the statistics of strong gravitational lensing of distant quasars by in- 
tervening galaxies, focusing on uncertainties in galaxy models (including velocity 
dispersion, luminosity functions, core radii and magnification bias effects) and 
on the parameters of the galaxy distribution and luminosity functions. In the 
process we derive new results on magnification biasing for galaxy lenses with non- 
zero core radii, and on how to infer the proper velocity dispersions appropriate 
for use in lensing statistics. We argue that the existing data do not disfavor a 
large cosmological constant. In fact, for a set of reasonable parameter choices, 
using the results of 5 optical quasar lensing surveys we find that a maximum 
likelihood analysis favors a value of fl D in the range ~ 0.25-0.55 in a flat universe. 
An open cosmology is not favored by the same statistical analysis. Systematic 
uncertainties are likely to be dominant, however, as these results are sensitive 
to uncertainties in our understanding of galaxy luminosity functions, and dark 
matter velocity dispersions, as well as the choice of lensing survey, and to a lesser 
extent the existence of core radii. Further observational work will be required 
before it is possible to definitively distinguish between cosmological models on 
the basis of gravitational lensing statistics. 

1 Introduction 

Recently there has been a renewed interest in the possibility that the cosmological constant 
may dominate the energy density of the universe, in order to resolve several cosmological 
conundra 0, |J, and more recently due to direct observational measurements of the distance 
redshift relation j|, For some time, one of the few cosmological pieces of evidence which 
has appeared apparently disfavoring this possibility has involved analyses of the statistics of 
gravitational lensing of quasars by galaxies || [7], ||, |9| . 

The effort to constrain cosmological parameters utilizing such statistics is of course highly 
dependent both on the quality of the existing data, and also on the robustness of the theo- 
retical inputs in the analysis. Because of the potential importance of the claimed constraints 
on a cosmological constant, it is worth re-examining in some depth the dependence of the 
resulting constraints on model dependent assumptions, as well as making some attempt to 
improve the models of galaxies one uses in the analysis. 
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In this paper, we re-analyze lensing statistics, concentrating on the role of the luminous 
E/SO galaxy distribution function and luminosity function parameters, and the self consistent 
modelling of galaxy parameters, including core radii. As we shall show, new considerations 
of the velocity dispersions of lensing galaxies || [H| , the core radii of lensing galaxies |7|, [TT 



and the magnification bias resulting from selection effects [g, £| [IT], [Tj] both affect the nature 
of the constraints one derives, and also imply that the constraints presently possible from the 
statistics of strong lensing may vary over a wide range. Nevertheless, we find for reasonable 
parameters, a best fit cosmology with a low matter density flat universe with a cosmological 
constant Moreover, we find that the data that favors a low Q Q mildly favors a flat 

cosmological constant dominated universe over an open universe. Five sets of quasar surveys 
are considered in order to compare predictions and observations. 



2 Theory layout: core radii, magnification bias, lensing 
statistics, and galaxy parameters 

In this section we first review the basic formalism we shall utilize in our statistical analyses, 
including the calculation of optical depths for lenses modelled as isothermal spheres with core 
radii. We then proceed to a new calculation of the relevant magnification bias which should 
be utilized in this situation. This is significant, because the trade off between reduction 
in lensing probability caused by finite cores is offset to some degree by an increase in the 
magnification bias due to lensing |J. The degree of this offset is important, however, if 
the two effects do not precisely cancel, as we find. We conclude this section with a short 
summary of our statistical maximum likelihood method. 



2.1 Optical depths and nonsingular isothermal spheres 

We begin by modelling the mass density distribution of elliptical galaxies with the following 
form pTf : 

where cxdm is the velocity dispersion of this system (presumably the dominant dark matter, 
a consideration we shall return to later), and r c is the core radius. When r c is zero, this 
model reverts to a singular isothermal sphere (SIS) model. With the density distribution 
from equation ([!]), one can calculate the lensing cross section, a cs , within which multiple 
images of gravitational lensing events are observed. 

a cs = Tralf(P) (2) 

where 

cap VotVts , , 

acT -R (l + z e )y 0S {6) 



3 As this paper was being finalized for submission a new investigation fl4|| also appeared which explores 
several similar issues and reaches compatible conclusions. 
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a cr is the critical radius, z is the redshift, and «o = 47r( £ ^ M ) 2 , is the bend angle for an 
isothermal sphere with core radius r c = 0. The Hubble constant H Q is 100 h km s _1 Mpc -1 . 
(Please note that this a cr was the size at lens redshift z\. It is not the critical radius at the 
present Universe.) The y quantities are angular diameter distances between the source, lens 
and observer, and will be discussed in more detail below. The quantity, j3 = r c /a CT as defined 
in Hinshaw & Krauss 01, and 



f((3) = 1 + 5/5- \f?- i^ + 4)i 



(4) 



These quantities arise in the determination of the bend angle, a, of the light trace from the 
source, which in the case of an isothermal sphere with a finite core radius, is a function of the 
velocity dispersion <tdm, impact parameter b, and core radius r c The general formula 
for the bend angle is: 



46 r°° 
a(b) = — / dr 

C Jb 



<9$ 



dr 



where 



$(r) 



d 3 r 



; Gp(r>) 



(5) 



(6) 



After calculating bend angles of multi images, one can then calculate the image separation, 
A9, and compare to measurements. 

The full lensing probability can then also be calculated as follows: 



T 



dr 



dz, 




FB{< m)f(P) 



Q (l + Zl ) 3 + Q R (1 + Zl ) 



(7) 



where 



c una 



o 



Ho 3 



The galaxy number density n is the luminous E/SO galaxy density, which is about 30 per 
cent of the total galaxy number density flS ], which is taken to be is 1.40 x 10~ 2 h 3 Mpc~ 3 , 
given by Loveday et al. [pl|. Spiral galaxies are not important in studies of strong lensing, 
because of their huge core radii 0. Qb. is the curvature term, and is the cosmological 
constant, so that Q Q + + Q\ = 1 in all cases. In a flat universe, yi s is simply y os -y i- 
However, in an open universe (with Q\ = 0): 



Vis = Vos 

y oi is the angular size distance fI7 |: 

1 



Vo 



sinh 



dz 



n (i + z) 3 + n R (i + zy + n A/ 



(9) 



(10) 



B{< m) is the magnification bias factor, which takes into account that lensed quasars are 
magnified, and thus have a larger probability of being observed than unlensed quasars, as we 
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shall discuss in the next section. We also have to integrate the lensing probability density 
over the lens luminosity distribution L, which is assumed to be given by the Schechter 
luminosity function 



^)^=(^)°exp(-^)rf(^) (11) 

The parameters of the Schechter function will be discussed in some detail later, as they 
induce perhaps the dominant uncertainties in the predictions for the gravitational lensing of 
quasars by galaxies. 



2.2 Magnification bias 



Magnification bias in a lensing analysis is due to the amplification of lensed quasar images. 
Because of this effect, an observer can observe a lensing event which cannot be seen if either 
it is not lensed or its magnification factor is not large enough to raise it above the minimum 
sensitivity of a flux-limited survey. A magnification bias factor, which enhances the optical 
depth compared to the bare optical depth calculated in the absence of this selection effect, 
can be estimated based on flux limits and the presumed luminosity distribution of quasars. 
This factor has been estimated earlier to be as large as 26 in the case of initial lensing 
surveys using a singular isothermal sphere (SIS) model for galaxy lenses |12|. Its formula 
and its value have been re-examined by Fukugita & Turner ||, and the value can be as low 
as 4, if the survey of a complete set of quasars with apparent magnitude less than 22 can be 
achieved. Amplification of a lensed quasar including a core radius of the lensing galaxy has 
been discussed in Hinshaw & Krauss 
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However, the bias factor has not been calculated 
before for this case. Here we carry out such an analysis. 

We begin with the definition of amplification (A) of a lensed quasar and also include the 
core radius in our formulas. Following Peebles [fTi 



.4 



Odd 




bdb 




xdx 


Oqd0Q 




Idl 




L\dL\ 



(12) 



where 9 is the angle between the lens and the pseudo image of the quasar, #q is the angle 
between the lens and the real quasar, b is the impact parameter, I is the transverse (projected) 
distance of the lens center from the line of sight, and x = b/a cr , L\ = £/a CT . We will see that 
the amplification A is a function of Li and (3 only, as discussed below. Following equation (6) 

< L x < L = 



in Hinshaw & Krauss |TT 
that equation to be 



-x 



^ < 1, and < < 
Vx 2 + (3 2 -f3 



X 



we can easily rewrite 



(13) 



This means that although A is a function of L±, x, and (3, we can rewrite A in terms of Li 
and (3 only. Using equation fll3]), we have 



dL x _ 1 
~dx~ ~ ^x 2 + (3 2 



Li 

— -2 

x 



(14) 
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and we can calculate the amplification A. For the total amplification, we sum up the absolute 
values of amplification caused by each individual image, i.e., 



cc ^ dec ^ 
L\dL\ 



(15) 



where Xi are the solutions of equation (|13"1), shown in Cheng & Krauss fl9|]. With 

1 



<^rr "I" 50 cr 7" c 



1 2 

— r 

2 c 



— r 

2 c 



1/2, 



+ 4a cr ) 3 / 2 



(we should remind our readers that L\ is actually identical to equation 
define the averaged amplification function of (3: 



<A>i 



Jo° dLi L\ 



— [ ° dLi Li A(Li, (3) 
Li Jo 



(16) 

above.) we can 
(17) 



The numerical values of < A > are listed in Table |l|. 

Our next step is to calculate the minimum value of amplification (Am) from the three 
quasar images for a given (3, which will be used in the magnification bias calculation. The 
numerical results are listed in Table 0. The values of /3, the corresponding minimum values 
of amplification, the corresponding L\ values at Am, and the corresponding L Q are listed in 
Table ^. We have also found that Am can be well fitted by 2/(L°- 65 ) and plotted in Figure |l[ 
We will use this approximation in our magnification bias calculation. 

We can now derive the magnification bias. We start with the probability density of 
amplification 

p(A)dA = 2A 2 M A~ 3 dA, A > A 



HI 



such that 



If we define 



roo 

/ p(A)dA = 1 



(19) 



A 



A r 



2.51og 10 A 
2.51og 10 A M 



then we have 



p{A)dA = p(A)dA = 0.8(ln 10)^l(T a8A dA 
The magnification bias is defined as: 

1 



B m) 



N Q (m) JA mi 



N Q (m + A)p(A)dA 



where 



. . / io a ( m " m °), if m < m Q 



(20) 
(21) 

(22) 

(23) 
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(Note the quantity 6 above should not be confused with the quantity used earlier to describe 
the impact parameter for lensing.) We then calculate the magnification bias: 

o|^y4|i 56 10( 6 - a ^ m - m °) if m a - A min <m<m Q 

ofe^M 6 if m Q < m 

(24) 

Interestingly, 6 has to be less than 0.8 from the above formula in order to have a reasonable 
B(m). The observational value of b is well below 0.8 and will be presented later. 

Next, as discussed in Fukugita & Turner 0, we have to average B{m) over the observed 
magnitude (m) distribution in order to obtain the relevant collective bias B(< m). We 
approximate a magnitude limited quasar survey using the selection function: 

„, \ _ f 1) when m is less (brighter) than the survey limit , 
1 0, when m is larger (dimmer) than the survey limit 

Ri , ^ - C, A " fiL„ dAp(A)N Q (m + A) 

B(< m) = jz ' ' 



Using equation (|23|) , we find 



a(a-fo) / 10 0.8(m-m o ) _ 10 0.8(m min -m o )\ gg 

S(< 77l) = ^ w ' -, r -, r- Ah Ah^ 

v ; (a -0.8) (0.8- b) \ 10 a ( m " m °) - 10 a ( m --- m °) J M a - 0.8 M 
when m min < m < m Q — A min < m 

B(< m) = ? — 1 ? [ + 10a(Wmfa>-mo)^5a 



I Q-^ Q 1 n fe(m-m ) ^2.56 Q ( a jO i n 0.8(m min -m o ) /|2 

+ 6(0.8-6) iU Am (a-0.8)(0.8-6) iU ^ 

when m m i n < m Q — A m j n < m < m 



0.8a / 10 fc ( m - m °) - l0 b ( m ^-m o ) \ ^ 



£(< m) = 7 — -, , -, r A 

v ' (o.8 — 6)6 y io a ( m - m °) — io a ( m n 1 m-m ) 

when m Q — A min < m min < m < m 

( 0.8 \ 10 fe (™-mo) _ lQb(m min -m ) 

B(< m) = -r— -, T7 ; Aif b 

v ' \ 8 _ \) b M _ lQa(m min -m )) i lQb(m-m ) _ I M 

a v ' 

when m Q — A min < m min <m <m 
B(<m) = - ? r-^ — ^ : (6-a + -^^10 a(mmi "- m ° ) ^ 



5 a 



M 



6(1 - 10 a ( m --- m °)) + a(10 fe ( m - m °) - 1) V a - 0.8 

I 0-^ a 1 Qb(m-mo) a 2.5b a6(a - 6) Q0,8(m min -m o ) 42 ] 

0.8-6 M (a-0.8)(0.8-6) M / 

when m min < m Q — A min < m Q < m (27) 
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2.3 Lensing statistics 



We utilize a maximum likelihood method based on Poisson statistics in our analysis [^0] 



If we have a complete and flux limited quasar survey, then we can calculate the expected 
number of lensing events, iV ex p, by summing up the probability of each quasar over all quasar 
samples, Nq, i.e., 



N, 



cxp 



i=l 



Zi) 



(28) 



Using equations (0) and fl28"|) , we can write down the expected number density distribution, 
Hi-, which is the integrand of equation (|7]). If we choose small bin sizes for source redshift, 
galaxy redshift, and luminosity jt, then there will be either no lensing events or one lensing 
event in each bin. Assuming appropriate Poisson probabilities Po,i and P\j, respectively, the 
likelihood function is 



C 



iVu 



0:t 



(29) 



where N \, s is the observed number of multiple-imaged lensing events, and N un is the number 
of the un-lensed events. After taking the logarithm of C, and taking the limit of summation 
to be an integral, we find the formula used in our lensing statistical analysis: 



ln£ = \nj\fj — iV ex p + constant 
i=i 



(30) 



with 
hi A/] 



= In 



In 



cPt 



dzi d (£) 



VolVls 



B(< m)f((3) 



,°"dm, 



fi G (i + Zl f + n R (i + Zl y + fi f 



u 

L* 



L* 



(31) 



+ constant 



3 Velocity dispersions and luminosity functions: dark 
matter and galaxies 

The probability of gravitational lensing can be seen, from equation (|7]) to depend strongly 
on both the luminosity function, and the velocity dispersion of galaxies. There is an explicit 
dependence on the latter in the lensing cross section. An implicit dependence arises as a 
result of the need to relate velocity dispersion to luminosity when performing the integral in 
equation (0). We have recently |21[] explored the model dependent considerations associated 



with extracting the relevant velocity dispersions of galaxies from the data, in the context of 
finite core isothermal models. We review the chief results here, along with discussing the 
interdependence of these estimates on the assumed luminosity distribution of galaxies. 

Returning to equation (|IJ), we should bear in mind that the velocity dispersion in this 
equation, which is assumed to be independent of radius, cannot be measured directly. What 
we can measure is the line-of-sight velocity dispersion o"i os (R), with the projected distance R 
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measured from the center of the observed galaxy (i.e., R is perpendicular to the line-of-sight 
to us). For a singular isothermal sphere, <tdm is not a function of R, but cr los (R) is, and this 
relationship is different when a galaxy has a finite core. Thus, one has to be careful how 
to use the measured velocity dispersion to derive the relevant quantity to utilize for lensing, 
namely whether it well approximates odm- Using a wide range of values of galaxy models, 
and self consistently solving the dynamical equations we demonstrated that of os /oD M 
is strongly sensitive to the core radius and galaxy anisotropy when the projected radius is 



less than 0.1i? e - This is consistent with observational data ||22|| . It can thus be dangerous 
to simply consider the line-of-sight velocity dispersion within 0.1R e as <tdm, even when core 
radii are small, as observations suggest they are for elliptical galaxies. However our results 
also indicated the following inequalities: 

1.16 < < 1-27 atR = 0AR e 

1.20 < < 1.30 at R = 0.54R e 

1.24 < < 1.37 atR = R e 

These indicate that the intrinsic scatter of <tdm will be less than 10 per cent if we can 
measure the line-of-sight velocity dispersion at R c or half R e . Then, we can multiply this 
velocity dispersion by the average value (e.g. 1.25 for R = 0.54R e ), in order to get <tdm- 
This argument is almost independent of the core radius of each galaxy and the anisotropy 
parameter f3. 

More important, we demonstrated that measuring the central velocity dispersion of a 
galaxy can give a misleading representation of <7dm> and in particular can give an overestimate 
of o"dm which can lead to a higher probability for lensing, and hence an inappropriately 
stringent bound on the cosmo logical constant. 

Finally, the above estimates are for the case of a purely finite isothermal distribution. 
If one adds to this distribution some central mass, such as a large central black hole, this 
will further increase the central velocity dispersion, and also change the relationship between 
<7dm and the line-of-sight velocity dispersion at R e and half R e . By measuring the velocity 
dispersion at both these points, however, one can hope to extract out the central mass 
contribution and also the isothermal contribution [F2TJ. (Alternately, it is clear that the 
velocity dispersion at R e will be less sensitive to the former contribution, and thus can be 
used to approximate the isothermal contribution.) The important point here is that the 
mean impact parameter for lensing at redshifts greater than unity is of order 1-10 kpc. At 
this range a central mass, if it is of the order of ~ 10 9 Mq, contributes for example at the 
5-10 percent level to the squared velocity dispersion, but it contributes merely at the 2-5 
percent level to the bend angle for lensing by the galaxy at this distance (the bend angle 
at some radius r is (2/n) times smaller for a point mass which has the same mass as an 
isothermal sphere at this radius, and hence produces the same velocity dispersion at this 
radius, because in the latter case the mass outside this radius contributes to the bend angle 
but not the velocity dispersion.) Since the fourth power of the velocity dispersion enters into 
the optical depth, interpreting a 10 per cent increase in the squared velocity dispersion using 
an isothermal sphere model predicts a 20 per cent increase in the optical depth, whereas the 
actual increase, if the contribution is from a central compact mass, is less than 10 per cent. 
Hence one must be concerned, when using velocity dispersions to estimate optical depths, 
about to what extent a possible central mass contributes to the former but not the latter. 
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To summarize these arguments one more time: The measured central velocity dispersion 
is not necessarily directly correlated to the quantity which is relevant for lensing statistics. 
This fact must be taken into account when attempting to utilize such statistics to constrain 
cosmological parameters. 

The last step in the determination of is the determination of the value of L* and 
the relation between L and odm- We consider the former first. First, to be consistent, it is 
important to calculate the luminosity function in the same band as one estimates velocity 
dispersions, and in which one performs lensing searches. Second, it is important to consider 
not the luminosity function for all galaxies, but rather for E/SO galaxies, as these dominate 
in the lensing statistics. An analysis of the latter has been carried out by Loveday et al. 
i"6fl , and yields, in the Bj band the relation, for the Schechter function: 



a = 0.20; M* = -19.71 + 51og 10 /i (32) 



In order to attempt to confirm this relationship in the Bt(0) band we calculated a lumi- 
nosity function utilizing data from a subset of sample in Faber et al. [^] and de Vaucouleurs 
et al. |24|]. In this case, we find a similar relation 



a = 0.15 ± 0.55; M^ t(0) = -19.66 + 5 log 10 h ± 0.30 (33) 

Note that this best fit M* differs from that used by Kochanek ||. At the same time we 
stress a fact emphasized by Kochanek, there is a complicated interplay in the determination 
of both M* and cTq M , so that one cannot arbitrarily vary one without the other. In any 
case, we utilize our relation above for the Bx(0) luminosity function in what follows below, 
although we also consider how our results would change if the parameters of the luminosity 
function change. 

In order to estimate the velocity dispersion of L* galaxies we found 42 galaxies with 
available velocity dispersions, 38 of which were appropriate for use in our analysis at R = R e , 



and 39 of which had reliable velocity measures at R = 0.54i? e When several different 
authors differed in their estimates of either R e or <7r c we utilized the weighted mean of the 
different estimates and incorporated the scatter in our error estimate. We then calculated the 



least-square- fit of the relation between log 10 (L/L*) and log 10 (<7R c /km s ) |2T[ and obtained 
log 10 (L/L*) = (-4.04 ± 0.49) + (1.89 ± 0.22) log^^/km s" 1 ) (34) 



If we set L = L* in equation fl3"4|), then we find the velocity dispersion of luminous 
elliptical galaxies at the effective radius, <7^ e = 135.9 ± 15 km s -1 . From the discussion 
above, we then multiply this number by 1.31, to get 0£, M ~ l-3l0R e ~ 178 km s -1 . 

Several issues are relevant to this result. In the first place, note that the L vs a relation in 



the above equation differs from the standard Faber- Jackson relation |[25|| . Note however that 
this FJ relation is appropriate for central velocity dispersions. In addition, it may be there 
case that the elliptical galaxies do not form a uniform population, but rather that bright 
galaxies and faint galaxies follow separate Faber- Jackson curves (this point was raised to us 
by J. Peebles, who has been investigating this issue |26|). It thus may not be appropriate to 
enforce this relation on the bulk sample. To explore this latter possibility, and because it is 
galaxies with the largest velocity dispersions which will dominate in the analysis of lensing 
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statistics, we considered a subset of our galaxy sample with log 10 (<7R e /km s x ) > 2.2 and 
rederived the L-a relation. In this case, we found: 

log 10 (L/L*) = (-6.20 ± 1.82) + (2.83 ± 0.79) log 10 (a fic /km s" 1 ) (35) 

In this case, we find a somewhat higher value of <7p M = 203 km s -1 as might be expected 
given the features of this subsample. It is worth noting that we have used the technique 
described earlier of using line-of-sight velocity dispersions at two values of R to extract out 
the asymptotic velocity dispersion which may be most appropriate for lensing, and find, 
coincidentally that for this subsample a value (T^ M = 179 km s _1 . Nevertheless, as a probe 
of the robustness of the gravitational lensing constraints we incorporate both of the above 
relations as well as the direct inferred values of (T^ M in our analysis. Note also that increasing 
the slope in the L vs <7dm relation has the same effect on lensing as decreasing a in the 
Schechter function, thus once again demonstrating the interdependence of these quantities 
in the deriving constraints. 

Finally, having derived two estimates for <jp M here, we must note that the lower value 
<Tq M pa 178 km s -1 is smaller than that utilized in several previous lensing statistical analy- 
ses. This may be due in part to what we claim is the inappropriate use of central velocity 
dispersions in other analyses, and while we reiterate that the value of a^ M cannot be inde- 
pendently varied in a self consistent analysis, and depends upon the form of the luminosity 
function, in order to allow for large possible systematic errors, we also used ct^m ~ 207 
km s _1 , which corresponds to a 2cr variation in the fit to <y^ M from first relation (42) above. 
Last, for purposes of comparison with Kochanek ||, we considered an even larger a^ M value, 
while at the same time adopting his choice of values for the Schechter function a — —1, L 
vs cr slope of 4, (and setting r c = in this case, see below) in order to explore the effects of 
variation in choice of M* and sample selection. 



4 Review of core radius 

The final quantity which we must determine in order to carry out our analysis is the rela- 
tionship between the core radius of E/SO galaxy and its luminosity. In Table |3], we list the 
inferred values of core radii from the Class I data in Lauer [?], distances of galaxies from 
Faber et al. and apparent magnitudes (Bt(0)) from de Vaucouleurs et al. P4| . With 
IVIbt(o) = -19.66 + 5 • log 10 h ± 0.30, log 10 (-^) for each galaxy in Table |3] can be calculated. 
We plot log 10 (-^-) against log 10 (r c /(/i _1 pc)) in Figure ||. With 13 samples in Table §, we 
have found the best fit is = (A) 1 ' 65 with r* = 45 h~ 1 pc and \ 2 = H- If we consider 
Ax 2 at the 95 per cent confidence level with three parameters of interest (i.e., r*,L*, and 
a power law relation between r c and L), such that the total x 2 = 18-8 (which is the 95 
per cent confidence contour line), with the power law relation and L* fixed, then we have 
r* = 45 ^Ijh^-pc. The best fit value of r* would in fact agree with an earlier investigation, 
if we neglect uncertainties from distances of galaxy samples, and normalize Mg , -) to the 
value chosen in Krauss & White |7J. 

An interesting test of this fit is to estimate the r c of NGC 7457, which is an SO galaxy 
in Lauer et al. P7[. Its apparent magnitude (Bt(0)) is 11.76 from de Vaucouleurs et 
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al. PU, and its distance is about 10 h 1 Mpc. Based on the formula in this section, we 



obtain r c = 5.2 h pc for the best fit of r*, and the 95 per cent confidence lower limit is 
r c = 3.2 /i pc. This result is marginally consistent with the lack of evidence for a core 
radius of NGC 7457 at a limit of « 3.4 pc (H G = 80 km s" 1 Mpc" 1 ) [|7|. 

In our analysis, with the parameters given above and later, without including magnifica- 
tion bias, we have found that the lensing probability is reduced due to the core radius by a 
factor 2 to 3, depending on the value of Q Q in flat cosmological models. With the magnifi- 
cation bias included, this reduction factor will be significantly decreased ||. However, one 
important result of our analysis is that these effects need not cancel out completely, depend- 
ing upon parameters of the luminosity function, so that as we shall see, the introduction of 
a finite core can in fact suppress the optical depth for lensing, even when magnification bias 
is taken into account. 



5 Quasar samples and lensing events 



5.1 Quasar samples 

We first utilized the Large Bright Quasar Survey (LBQS) catalogue of Hewett, Foltz & Chaffe 
||28|| . There are 1055 quasar samples in that catalogue. All quasars are between redshift 0.2 
and 3.4, and their apparent magnitudes (Bj) are between 16 and 18.9. We include these flux 
limits to calculate the magnification bias in our statistics. There was only one lensing event, 
LBQS1009-0252, observed in this survey. 

As a comparison, we also used the Snapshot survey by Maoz et al. 



29 combined with 



other surveys by Crampton, McClure & Fletcher |30[| , Surdej et al. fTTH , and Yee, Filippenko 
& Tang ||32|| . There are 648 quasars in this total sample set. There were 4 lensing events 
in these surveys. In modelling the magnification bias for this set we utilized the flux limit 
of 19.5, appropriate for the Snapshot survey, which contains the largest number of quasars, 
slightly more than the Surdej et al. survey. Our results do not depend sensitively on this 
choice, however, and in any case the flux limit in the latter survey is similar to that in the 
Snapshot survey. 

Finally, if we consider all five surveys together, then we will have 1615 quasars in total, 
with five lensing events. In our analysis of this combined sample we again utilized the flux 
limit from the Snapshot survey (lower than that for the LBQS survey) to calculate the 
magnification bias. 



5.2 Lensing events 

Kochanek []] has previously discussed the rationale for considering the 5 strong gravi- 
tational lensing events mentioned above among all optical gravitational lensing candidates 
for use in the calculation of lensing statistics. We list the quasar redshifts, lens redshifts, 
and angular splittings for these events in Table f§ |33], R41 EB, |36], R7L [M B9|. In order to 



apply our statistical model, we actually need to know the redshift of each lens. At present, 
only the redshifts of lenses in two systems (Q0142-100 and PG1115+080) have been firmly 
measured. For those lenses which lack galaxy redshift information, we chose the most likely 
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absorption lines as lens redshifts. For systems H1413+117 and Q1208+1011, we used two 
possible redshifts in our analysis, as two different lenses are possible candidates. 

With the information in Table [|, we can calculate the likelihood function in equation (|32|). 
What we need are the core radius and velocity dispersion for each galaxy. Recall the angular 
separation is well approximated by A9 = 2^1 — 2(3 cto^ [0- Using the results of section 
3, we employ power law relations between luminosity and velocity dispersion as either = 
(£em)7 ; wit h (i) 7 = 1.89, cr* = 178 km, s" 1 or (2) 7 = 2.83, a^ M = 203 km s _1 . This 

"dm 

implies a power law relation between core radius and velocity dispersion, i.e., = (£em) 3 - 12 

r c ff DM 

or % = (2em) 4 -67_ Combining this equation with a known A8 value, and with a given 

r c "DM 

cosmological model, in principle, we can solve r c and <tdm for each lensing galaxy. Actually, 
except for the lensing galaxy for the system Q1208+1011, all other four lensing systems 
allow a second solution with <tdm at least 1400 km s -1 , which is not a suitable scale for a 
galaxy. In our analysis we select the smaller but reasonable velocity dispersion. In Tables |5] 
and §, we list our theoretical predictions of odm, ^c, and B T (0) for the lensing galaxy for 
each lensing system in different cosmological models. Note that the inferred o"dm values are 
reasonable, but have an average which exceeds the value of a^ u we have used. This is to 
be expected. Lensing galaxies will have preferentially larger velocity dispersions than the 
mean, simply because these galaxies are weighted more heavily in the probability function. 
It would thus be inappropriate to use the mean for lensing galaxies to model all galaxies. 
For the five events considered here, we find that that (or j^) increases when VL in- 

"dm l 

creases for all five events in flat universe models. This can be explained as follows: y os j (y iyis) 
increases when Q Q increases for a fixed source redshift and a fixed lens redshift [[?], 40 . For 
a fixed (observed) A9, when y os /(yoiyis) increases, has to increase in order to prevent 

"dm 

(3 from approaching 1/2. The ratio as a function of VL seems to increase faster when 

°"dm 

AO is smaller. In fact, there is no solution for r c and <tdm of the system Q1208+1011 when 

Q D > 0.95, if Zi = 2.9157 is chosen. This is because — decreases when Q increases for a 

yos 

fixed source redshift and a fixed lens redshift. As a result, A6 from the system Q1208+1011 
is too 'large' in the Q Q = 1 universe model if the given source redshift and lens redshift are 
accurate and correct. If more 'large' angular splitting lensing systems are observed at high 
redshifts, this will more strongly favor a low matter density universe, although of course the 
number of lensing events will play an important role in constraining models. In open universe 
models, the results are similar to flat universe models, but the £ P M ratio as a function of fi Q 

""dm 

is smoother. 



6 Numerical results 
6.1 Magnification bias 

With the parameters discussed above, we can calculate the lensing optical depth. One 
interesting task is to compare our bias factor to the SIS model in Fukugita & Turner 0. We 
define the average of B{< m) as following: 

<B(< m)> = - / B(< m)dr (36) 
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where r is the optical depth of the lensing events. In addition to the parameters discussed 
previously, We choose m min to be 16, m Q = 19.15, a = 0.86, b = 0.28 from Hartwick & Schade 
PH . In Figure || we have plotted <B(< rri)> as a function of quasar (source) redshift in 
different cosmological models, with m = 18.9, 19.5, and 22. 

Although we set a lower limit for the integrals in equation (|26"D, there is no significant 
impact if we choose m m j n to be either 16 or smaller (such as negative infinity as in Fukugita 

6 Turner ||). In Figure |3|, we have shown that <B(< ni)> is almost a constant at higher 
source redshift (z s > 1) in each cosmological model. Although the optical depth of the 
Q\ = 1 flat universe model is larger than the optical depth of the Q Q = 1 universe model, the 
averaged total magnification bias goes in the other direction. This is also true in the open 
universe model. Comparing to the bias values in Fukugita & Turner |J, we have <B{< m)> 
around 13 when m = 18.9, but they have 7.33. We have <B(< m)> around 7, but they 
have B(< m) = 4.25, when m = 19.5. We have <B(< m) > around 4, and they have 
B{< m) = 2.63, when m = 22. These numbers confirm that although adding a core radius 
in the gravitational lensing analysis will generally reduce the lensing probability by a factor 
2 to 3, the magnification bias including the core radius is about 1.6 times (depending on the 
flux limit of the quasar survey) larger than the SIS model. This is due to the fact that with 
a non zero core radius the minimum value of the total amplification of all quasar images is 
generally larger than 2, the value for the SIS model. 

6.2 Flat universe models 

Incorporating the magnification bias results discussed above we first consider the statistical 
analysis in which quasars from all 5 surveys were utilized. For flat universe models, these 
are presented in columns 2 to 7 in Table [F[ As mentioned earlier, the apparent magnitude of 
the combined quasar survey limit has been set to be 19.5. If we consider our best fit galaxy 
parameters it is clear from columns 2 to 5 that the expected number of lenses is a good fit, 
and the likelihood function is minimized, when Q Q « 0.2. The strict 95 per cent confidence 
level (for 1 degree of freedom) on this quantity for these galaxy parameters ranges from 
0.07 < Q Q < 0.55 (column 3) or 0.07 < Q Q < 0.65 (column5). However, if we consider the 
intrinsic scatter merely in velocity dispersions and core radii, then we find a different result. 
For example, consider a smaller r*, 28 h~ 1 pc, and a larger a*, 207 km s _1 . Columns 6 and 

7 show a best fit around Q D = 0.5, and the 95 per cent confidence level is Q D > 0.24. This 
result here is comparable, although less stringent than the result (with Q Q = 1 as the best 
fit value) given in Kochanek ||, without including core radii effects. This demonstrates the 
sensitivity of existing statistical constraints to the values of the assumed velocity dispersions 
and core radii. Note that the inclusion of core radii, and the difference in the assumed value 
of M* presumably account for the reduction in the minimum allowed value of fl D compared 
to Kochanek 0. Also note that inclusion of core radii is largely responsible for shifting the 
peak of the likelihood function in this case toward values lower than Q — 1. In any case, 
while suggesting that Q < 1 is preferred by our best fit parameters, the divergence in results 
demonstrates the sensitivity of existing statistical constraints to the values of the assumed 
velocity dispersions and core radii. The likelihood functions mentioned above are plotted in 
Figure f|. 

Setting Q = 0.2 in Figure we plot the expected number of lensing events and the 
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observed (five) events vs source redshift (z s ) . The shape of the predicted curve is in agree- 
ment with the observation, given the very limited statistics. If one had a complete quasar 
survey up to a high redshift (larger than 2.5), one could use this shape function to probe 
more clearly the agreement, or lack thereof, between theory and observation. 

In columns 8 to 13 in Table [7|, we list results obtained by using only the LBQS survey and 
the single lens candidate, LBQS 1009-0252. Because only one lens is considered, the results 
more strongly favor a higher Q universe. Nevertheless, while the best fit of fl is around 
0.8, the 95 per cent confidence level lower limit is about Q > 0.15. The large width of the 
allowed region presumably reflects both the reduced statistics, and the effect of non-zero 
core radii. The likelihood functions are plotted in Figure |6|. Columns 12 and 13 involve the 
parameters used in Kochanek || (i.e., a = —1, 7 = 4, r* = 0, a* = 225 km s _1 ). The best 
fit of fl from the likelihood analysis is 0.95, with 95 per cent confidence level, Q > 0.2 in 
flat universe models (see also Figure ||). This result is consistent with the result given by 
Kochanek |J, but again the lower limit on Q Q is reduced. 

For comparison, in columns 14 to 16 the result using the four combined quasar surveys 
PO] , |3~T| , |32j without LBQS are shown. As expected this tends to favor slightly smaller Q Q 



than the five survey analysis. The best fit of Q Q is about 0.15. The 95 per cent confidence 
levels are 0.03 < Q Q < 0.43 and 0.03 < Q Q < 0.38 for the last two columns. Statistical values 
are not available at high fl values in the last two columns, because <tdm arid r c of the lensing 
galaxy in system Q1208+1011 cannot be fitted, as discussed in the previous section. 

As another probe of the galaxy lensing parameters, if we define the mean angular splitting 

as 

<A9>=-f A9dr (37) 



with A6 = 2 v /T^2/3 |9[, we find that < A6 >^ 279 (7 = 1.89, o£ M = 178 km 



s- 1 ) or 2"0 (7 = 2.83 

,fDM = 203 km s when the source redshift is less than 5 in flat 

universe models. However, if we calculate the mean value of the angular separation from 

Table f|, we find a mean value l'/542. From the formula of < Ad > for the SIS model 

||, one can reduce the mean angular splitting by either decreasing o* or increasing the 

value of the parameter 7 (for fixed Schechter a). Indeed, if we choose cx^m = 178 km s -1 

and 7 = 2.83, the combination suggested by extracting away the central mass contribution 

from the asymptotic velocity dispersion as discussed previously, we find < A9 >k, l'/5. From 

Figure |], one can infer in this case a best fit for Q Q around 0.23, with a 95 per cent confidence 

interval around 0.07 < Q Q < 0.6. 

6.3 Open universe models 

We repeat the above analysis for open universe models. When all 5 surveys are included in 
the analysis the results are displayed in columns 2 to 5 in Table [|. In this case the best fit 
is fl ~ 0. Note, that the expected number of lensing events is smaller than the observed 
number when one considers more realistic values Q > in this case. The 95 per cent 
confidence limits are nevertheless fairly broad, with either Q Q < 0.65 (column 3) or Q < 1 
(column 5). 



We next consider just the four quasar surveys (Crampton et al. P0f , Maoz et al. [^9 



Surdej et al. [plj| , and Yee et al. |j32[). The results are listed in columns 6, 7, and 8 in 
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Table The best fit value is again Q Q = 0, with a 95 percent confidence limit of Q Q < 0.3 
in this case. 

Finally, if we consider the LBQS survey alone, with a single lensed event, it is clear that 
higher density open universe models will fare better. Since in fact the best fit value of Q a is 
in this case close to 0.95, the open universe model likelihood function in this region is then 
very similar to that for the flat model. As a result, the likelihood function and predicted 
number of events will have the same behavior (near Q = 1) as the values shown in columns 
8-13 of Table [?| and displayed in Figure |] for the flat universe case. As a result, we do not 
explicitly re-display the values again in Table || 

In general open universe models tend to predict a redshift distribution for the lenses 
which is in poorer agreement with observation, although again because of the very limited 
statistics this is not a quantitative problem at this point. It is also worth noting that one 
can increase the number of lensing events by increasing <J^ M as seen in column 4, but in 
this case the mean predicted angular splitting will be increased, as we have discussed in the 
previous subsection. 

7 Conclusions 

The various different fits described here in Tables [7] and |8] demonstrate that the constraints 
one derives on cosmological models from existing lensing statistics depend in detail on the pa- 
rameters one uses to model galaxy distributions. Nevertheless, the model parameters which 
we suggest are favored involve a low density, cosmological constant dominated universe — the 
model which coincidentally is favored by other astrophysical data at the present time. In 
our analysis, we have been careful to self consistently obtain the various lensing parameters 
(velocity dispersions, core radius, magnification bias, etc.). As a result, the general features 
of our analysis are expected to be robust, even though the intrinsic scatter in values of the 
selection functions, galaxy velocity dispersions and galaxy core radii may be large at the 
present time. Although the selection function in this paper is as simple as in Fukugita & 
Turner ||, a real selection function should at least account for the image separation and 
flux ratio. If we include these factors in our paper, they will reduce the number of pre- 
dicted lensing events. With the observed number of lensing events remaining unchanged, 
this means that the lensing statistics will favor a flat universe with a larger cosmological con- 
stant value — in even greater disagreement with Kochanek's conclusion |§. Quantitatively 
our analysis suggests a best fit value of Q Q ranges between 0.25-0.55 in a flat universe if five 
lensing events in the 5 optical surveys considered here are incorporated. (Note however that 
this value can raise to 0.95 if only one lensing event and a single survey are included in the 
analysis.) The distribution of predicted events is in a good agreement with the distribution 
of observed events in source redshift space (see Figure |5]). Note however that the results 
depend on which lensing surveys are utilized and thus the number of lenses used. Clearly 
what is required in order to distinguish between model parameters are greater statistics from 
a complete full sky quasar survey up to high redshift. 

The predicted number of lenses in open universe models do not vary as sharply as flat 
universe models do over the range of Q D , so while the preferred value of Q Q is around 
if lensing events from all five lensing surveys considered here are incorporated, the allowed 
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range remains very broad. The fact that open universe models tend to predict too few events 
when compared to 5 lensing events in the five surveys has also been demonstrated by Chiba 



& Yoshii [13]. Note that even if a larger o"5 M is chosen compared to those utilized here, in 
order to increase the number of lenses, then the mean angular splitting will become larger 
than the mean value of the observations in this case. Finally, open universe models tend 
to predict more lower redshift lenses than flat universe models, and this may be a useful 
discriminant in the future. 

We remind readers that the lensing optical depth is proportional to the fourth power of 
the velocity dispersion. A 10 per cent change of velocity dispersion can cause 46 per cent 
change in the optical depth. Thus, better estimates of velocity dispersions for E/SO galaxies 
are necessary before the statistical limits can improve. We have argued here that determining 
the appropriate velocity dispersion for use in lensing statistics, when core radii are not zero 
can be somewhat subtle. Self consistent modelling of the gravitational potentials may be 
necessary before line-of-sight velocity dispersions can be used to infer dark matter velocity 
dispersions. In addition, the effect of compact central mass conglomerations may bias the 
interpretation on extracts from using central velocity dispersions in the derivation of lensing 
statistics. 

The parameter 7, relating velocity dispersion and luminosity is also important in grav- 
itational lensing studies. A larger value of 7 reduces the mean angular splitting and also 
pushes Q Q toward a low value and a narrow 95 per cent confidence level. (These effects can 
be easily understood from equations (|28| ) and (|3"2D.) Therefore, improved observations of 
E/SO galaxy velocity dispersions along with luminosity measurements in order to reduce the 
intrinsic scatter in the data as well as to better model the Schechter parameters would be 
useful. 

We have also demonstrated that adding a core radius in the lensing analysis will generally 
reduce the total optical depth, even though the magnification bias factor will be higher than 
that predicted SIS models. The overall reduction of the lensing optical depth for non zero 
core radii is parameter dependent, however. For example, if we choose a steeper power 
law relation between the velocity dispersion and the luminosity (L oc cr 4 ), with Schechter 
a = — 1, then including a core radius in the analysis produces no significant change in the 
number of expected lensing events, and produces shifts in the best fit likelihood analysis for 
f2 Q of < 0.1. With a shallower power law relation and the Schechter parameters used in this 
paper, we have however found a 30 per cent to 40 per cent change of the number of predicted 
lensing events by adding a core radius. This result is easily understood: a shallow power law 
effectively more strongly weights larger core radii when L is less than L* . 

Although the apparent magnitudes are given in Bt(0) band without K-corrections in 
Tables |5] and |6], these values nevertheless suggest why observers have not yet found the 
lensing galaxy in systems H1413+117, LBQS1009-0252, and Q1208+1011. We also suggest 
examining system Q1208+1011 carefully, because this event, because of its angular splitting 
and assumed the high lens redshift, helps drive the fits to a low Q universe. Note that 
recent studies of the faint lens image of H1413+117 suggest a value of the lens brightness 



consistent with that calculated in Table [| ] p2| , 

To summarize: the self consistent incorporation of core radii in galaxy models, and fitting 
the distribution of galaxies to the E/SO galaxies which dominate lensing statistics combine 
together to suggest a statistical best fit of predictions to the results of existing optical quasar 



16 



lensing surveys for a flat universe model with Q Q in the range 0.25-0.55. Considerable system- 
atic uncertainty persists however in the appropriate velocity dispersion to use in models, as 
well as the luminosity- velocity dispersion relation, and the determination of the appropriate 
luminosity function parameters. In addition, how to appropriately combine the results of 
different surveys by different groups is not obvious. Since the results depend upon all these 
factors, it may be premature to argue definitively in favor of this best-fit scenario, although 
it is encouraging that it agrees with that obtained using other cosmological observables. In 
any case, we have derived here several new results here associated with magnification biasing 
and galaxy models with core radii, as well as the determination of the appropriate velocity 
dispersion to use in galaxy lensing models. As data improves, we expect that the techniques 
described here will be useful in further constraining theoretical models. 
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Figure 1: Plot of natural logarithm of A M versus parameter (3. The dots are numerical 
results, and the solid curve is the best fit curve: 2/(L°- 65 ). 
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Figure 2: The best fit plot of log 10 jj vs log 10 (r c /(/i _1 pc)). There are 13 samples with error 
bars shown in the figure. The middle solid line is the best fit, and the two dash lines give 
the 68 per cent confidence level, assuming the best fit power law. The two dash-dot lines 
give the 95 per cent confidence level under the same condition. 
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Figure 3: Plots of <B(< m)> versus quasar (source) redshift. The solid curve sets are for a 
flat universe model, with Q Q = 1. The dash-dot curve sets are also for a flat universe model, 
but with Q = 0, and Q\ = 1- The dot-dot curve sets are for an open universe model, 
with Q Q = and zero cosmological constant. The black (highest) curves represent quasar 
survey limit, m = 18.9. The gray curves represent quasar survey limit, m = 19.5. The light 
gray (lowest) curves represent quasar survey limit, m = 22. We choose a^, M = 178 km s _1 , 
r* = 45/i _1 pc, and 7 = 1.89 to generate these curves. 
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Figure 4: Plots of likelihood functions vs Q Q . The solid curve is the plot of the third column 
in Table 0, with r* = 45 /i _1 pc, <j^, m = 203 km s^ 1 , and 7 = 2.83. The dash-dash curve 
represents the result for r* = 45 /i _1 pc, <j^ M = 178 km s -1 , and 7 = 1.89. The dot-dot curve 
is the plot for r* = 28 h^pc, = 207 km s -1 , and 7 = 1.89. 
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Figure 5: The predicted number of events and observed events are plotted as a function 
of quasar redshifts from all 5 quasar surveys with fl Q = 0.2 in a flat universe model. The 
small dots making up the curve give the total predicted number of lensing events within a 
given source redshift. The ellipses give the number of observed events. The parameters are 
cr* M = 203 km s" 1 , r* = 45k- 1 pc, and 7 = 2.83. 
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Figure 6: Plots of likelihood functions vs Q Q . The solid curve is the plot for r* = 45 h _1 pc, 
cr^M — 203 km s _1 , and 7 = 2.83 in flat universe models using only the LBQS result. The 
dash-dash curve represents the plot for r* = 45 ft, -1 pc, a^ M = 178 km s -1 , and 7 = 1.89. 
The dot-dot curve is the plot for r* = 0, a^ M = 203 km s -1 , 7 = 4 and luminosity function 
parameter a = — 1. 
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(3 0.025 0.05 0.075 0.1 0.125 0.15 0.175 0.2 0.225 
<A> 4 5.81 7.08 8.44 9.99 11.9 14.1 16.9 20.4 24.9 



(3 0.25 0.275 0.3 0.325 0.35 0.375 0.4 0.425 0.5 
<A> 30.7 38.7 49.9 66.3 91.6 134 213 383 oo 



Table 1: Averaged amplification as a function of (3. 
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0.137874 


0.168375 




0. 


3 


32.6748 


0.0910677 


0.114322 




0, 


35 


60.7457 


0.0549151 


0.0708916 




0. 


4 


141.863 


0.0278588 


0.037016 




0. 


45 


585.195 


0.00920176 


0.0126015 




0. 


49 


14931.9 


0.000780135 


0.001096 




0. 


5 


oo 








Table 2: 


A M , Li(A m ), and L Q as a 


function of (3. 


NGC 




e c 


(arcsec) 


Distance (/i _1 Mpc) m(B T (0)) 


720 




4.58±0.01 


20.5±4.4 


11.13 


741 




1.91±0.03 


53±11 


12.18 


1395 




U 


36±0.12 


19.9±1.9 


10.52 


1407 




3.21±0.07 


19.9±1.9 


10.51 


1600 




3.73±0.09 


40.2±4.9 


11.85 


3379 




1.66±0.07 


8.6±1.3 


10.17 


4261 




2.61±0.05 


27.8±5.9 


11.32 


4365 




2.33±0.02 


13.33±0.71 


10.42 


4374 




2.99±0.02 


13.33±0.71 


9.91 


4472 




3.7±0.03 


13.33±0.71 


9.26 


4636 




3.24±0.17 


13.33±0.71 


10.37 


4649 




3.41±0.11 


13.33±0.71 


9.7 


5846 




2.51±0.14 


23.4±2.8 


10.87 



Table 3: Core radii of sample galaxies and their corresponding distances and apparent mag- 
nitudes. 
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Source 


z s 


Zl 


mi 


A9 (arcsec) 


Q0142-100 


2.719 


0.49 


R « 19 


2"2 


PG1115+080 


1.722 


0.294 


R « 19.8 


2"3 


H1413+117 


2.551 


1.4382" (1.6603 a ) 




l'/23 


LBQS1009-0252 


2.739 


0.869 a 


R > 21 


l'/53 


Q1208+1011 


3.803 


1.1349 a (2.9157 a ) 




0"45 



a The assumed galaxy redshifts, based on the most possible absorption 
lines by observers in literature. 



Table 4: Parameters of gravitational lensing systems. First column: the name of each quasar 
in each lensing system. Second column: the source redshift. Third column: the lens redshift. 
Fourth column: the apparent magnitude (in R band) of each galaxy, if the value is available. 
Fifth column: the angular splitting of each lensing system. We take the average separation 
value of the quadruple system H1413+117. 







tt = 


=0, tt A =l 






=0.2, Q A =0.8 




source 


<7 DM (km s 




r c (h- l pc) 


B T (0) 


<7 DM (km s" 




B T (0) 


Q0142-100 


217 




62 


21.8 


233 


86 


21.4 


PG1115+080 


221 




67 


20.4 


230 


81 


20.1 


H1413+117 


223 




69 


25.2 


275 


185 


23.9 


H1413+117 a 


250 




119 


25.3 


324 


398 


23.8 


LBQS1009-0252 


198 




40 


23.9 


224 


72 


23.1 


Q1208+1011 


106 




2.1 


26.7 


124 


4.5 


25.7 



Galaxy redshift is 1.66. 



Table 5: Theoretical predictions of velocity dispersions, core radii, and apparent magnitudes 
of lenses in flat universe models. K-correction is not included in the fourth and seventh 
columns. We choose a^ M = 203 km s -1 , r* = 45/i _1 pc, and 7 = 2.83. 
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n Q = 


=o, n R =i 




n D -- 


=0.25, tt R =0.75 




source 


<tdm (km s _i ) 


r c (/i _1 pc) 


B T (0) 


o- DM (km s" 


_1 ) r c (/i _1 pc) 


B T (0) 


Q0142-100 


254 


129 


21.0 


254 


127 


20.9 


PG1115+080 


244 


107 


19.8 


244 


107 


19.8 


H1413+117 


316 


354 


23.3 


320 


376 


23.1 


H1413+117 a 


390 


954 


23.1 


420 


1338 


22.7 


LBQS1009-0252 


252 


123 


22.5 


251 


122 


22.4 


Q1208+1011 


141 


8.3 


25.1 


140 


8.0 


24.9 



a Galaxy redshift is 1.66. 



Table 6: Theoretical predictions of velocity dispersions, core radii, and apparent magnitudes 
of lenses in open universe models. K-correction is not included in the fourth and seventh 
columns. We choose <r DM = 203 km s" 1 , r* = 45/i _1 pc, and 7 = 2.83. 



25 



cxp 



a, 6 



Kf a 4 

cxp 



cxp 



cxp 



In £ c '3 In £ c '3' /1 



ln£ 



i\T a ' c 

exp 



ln£ a ' 



ln£ 



ln£ 6 



N c ' e 

cxp 



ln£ c 



cxp 



to 





0.05 

0.1 

0.15 

0.2 

0.25 

0.3 

0.35 

0.4 

0.45 

0.5 

0.55 

0.6 

0.65 

0.7 

0.75 

0.8 

0.85 

0.9 

0.95 

1. 



21.0 

10.8 

7.5 

5.8 

4.7 

3.9 

3.4 

3.0 

2.6 

2.4 

2.1 

2.0 

1.8 

1.7 

1.5 

1.4 

1.3 

1.2 

1.2 

1.1 

1.0 



-28.16 
-19.93 
-18.11 
-17.56 
-17.47 
-17.59 
-17.74 
-18.03 
-18.36 
-18.70 
-19.04 
-19.39 
-19.74 
-20.08 
-20.42 
-20.76 
-21.09 
-21.42 
-21.74 
-22.05 
-22.36 



22.2 

11.3 

7.9 

6.1 

4.9 

4.1 

3.6 

3.1 

2.8 

2.5 

2.2 

2.0 

1.9 

1.7 

1.6 

1.5 

1.4 

1.3 

1.2 

1.1 

1.1 



-27.95 
-18.98 
-16.86 
-16.11 
-15.87 
-15.80 
-15.93 
-16.12 
-16.36 
-16.61 
-16.88 
-17.16 
-17.43 
-17.70 
-17.98 
-18.24 
-18.51 
-18.76 
-19.02 
-19.26 
-19.51 



47.7 

24.8 

17.4 

13.5 

11.0 

9.3 

8.0 

7.1 

6.3 

5.7 

5.1 

4.7 

4.3 

4.0 

3.7 

3.5 

3.2 

3.0 

2.9 

2.7 

2.6 



-54.69 
-33.43 
-27.25 
-24.30 
-22.63 
-21.63 
-20.99 
-20.59 
-20.34 
-20.10 
-20.03 
-20.02 
-20.04 
-20.09 
-20.16 
-20.25 
-20.35 
-20.46 
-20.58 
-20.71 
-20.84 



14.7 

8.3 

6.0 

4.7 

3.9 

3.3 

2.9 

2.5 

2.3 

2.0 

1.8 

1.7 

1.6 

1.4 

1.3 

1.3 

1.2 

1.1 

1.0 

1.0 

0.9 



-14.91 
-8.95 
-6.92 
-5.87 
-5.24 
-4.83 
-4.43 
-4.24 
-4.10 
-4.01 
-3.94 
-3.89 
-3.86 
-3.84 
-3.83 
-3.83 
-3.83 
-3.84 
-3.85 
-3.87 
-3.89 



15.5 

8.8 

6.3 

5.0 

4.1 

3.5 

3.0 

2.6 

2.4 

2.1 

1.9 

1.8 

1.6 

1.5 

1.4 

1.3 

1.2 

1.1 

1.1 

1.0 

1.0 



-15.49 
-9.13 
-6.95 
-5.81 
-5.13 
-4.56 
-4.25 
-4.03 
-3.87 
-3.76 
-3.67 
-3.61 
-3.56 
-3.53 
-3.51 
-3.50 
-3.49 
-3.49 
-3.49 
-3.50 
-3.51 



19.5 

11.3 

8.3 

6.6 

5.4 

4.7 

4.1 

3.6 

3.2 

2.9 

2.7 

2.5 

2.3 

2.1 

2.0 

1.8 

1.7 

1.6 

1.5 

1.5 

1.4 



-19.16 
-11.53 
-8.88 
-7.50 
-6.65 
-6.09 
-5.70 
-5.41 
-5.20 
-5.05 
-4.93 
-4.84 
-4.77 
-4.72 
-4.68 
-4.66 
-4.64 
-4.63 
-4.62 
-4.62 
-4.62 



16.1 

7.7 

5.2 

4.0 

3.2 

2.7 

2.3 

2.0 

1.7 

1.6 

1.4 

1.3 

1.2 

1.1 

1.0 

0.9 

0.9 

0.8 

0.7 

0.7 

0.7 



-20.63 
-15.37 
-14.57 
-14.55 
-14.78 
-15.04 
-15.43 
-15.84 
-16.25 
-16.67 
-17.08 
-17.50 
-17.92 
-18.35 
-18.78 
-19.22 
-19.70 
-20.22 
-20.85 
N/A 
N/A 



-20.73 
-15.75 
-15.12 
-15.20 
-15.52 
-15.85 
-16.29 
-16.75 
-17.22 
-17.68 
-18.14 
-18.60 
-19.06 
-19.52 
-19.99 
-20.48 
-20.99 
-21.54 
-22.21 
N/A 
N/A 



a All five quasar surveys and lensing events are used. 
b ct* m = 203 km s _1 , r* = 45ft- 1 pc, and 7 = 2.83. 
c a* M = 178 km s _1 , r* = 45/i" 1 pc, and 7 = 1.89. 
d a ^ M = 207 km s _1 , r* = 28/i" 1 pc, and 7 = 1.89. 
e Only LBQS and LBQS 1009-0252 are considered. 
f ct^m = 225 km s -1 , r* = 0, 7 = 4, and a = —1. 



9 Based on surveys by Crampton et al. [30], Maoz et al. [29], Surdej et al. [31], and Yee et 
al. [32]. Four lensing events are used. z t = 2.9157 in system Q1208+1011 is used. There is 
no solution for r c and ct D m of the system Q1208+1011 when fi D > 0.95. 

h zi = 1.66 in system H1413+117. 

Table 7: Expected lensing events and likelihood analysis. This table is for flat universe models, i.e., fi Q + Qa = 1. A^ X p is our 
theoretical prediction of the number of lensing events, and ln£ is the logarithm of the maximum likelihood function. z\ = 1.438 
in system H1413+117 and z\ = 1.1349 in system Q1208+1011 are generally used, unless otherwise noted. 





pra,b 
cxp 


\nC a ' b 


/ya,c 
cxp 


In C a ' c 


N d 

cxp 


ln£ d 


In C d ' e 





2.1 


-19.43 


5.0 


-22.34 


1.3 


-16.04 


-16.99 


0.05 


1.9 


-19.61 


4.7 


-22.38 


1.3 


-16.40 


-17.39 


0.1 


1.9 


-19.79 


4.5 


-22.43 


1.2 


-16.74 


-17.76 


0.15 


1.8 


-19.96 


4.3 


-22.49 


1.2 


-17.06 


-18.12 


0.2 


1.7 


-20.06 


4.1 


-22.49 


1.1 


-17.29 


-18.38 


0.25 


1.6 


-20.23 


4.0 


-22.56 


1.1 


-17.59 


-18.70 


0.3 


1.6 


-20.39 


3.8 


-22.64 


1.0 


-17.87 


-19.01 


0.35 


1.5 


-20.55 


3.7 


-22.72 


1.0 


-18.14 


-19.31 


0.4 


1.5 


-20.70 


3.6 


-22.80 


0.9 


-18.41 


-19.61 


0.45 


1.4 


-20.86 


3.5 


-22.89 


0.9 


-18.68 


-19.89 


0.5 


1.4 


-21.01 


3.3 


-22.98 


0.9 


-18.94 


-20.18 


0.55 


1.3 


-21.15 


3.2 


-23.07 


0.9 


-19.21 


-20.47 


0.6 


1.3 


-21.30 


3.2 


-23.16 


0.8 


-19.48 


-20.76 


0.65 


1.3 


-21.44 


3.1 


-23.25 


0.8 


-19.75 


-21.06 


0.7 


1.2 


-21.58 


3.0 


-23.34 


0.8 


-20.04 


-21.36 


0.75 


1.2 


-21.71 


2.9 


-23.43 


0.8 


-20.35 


-21.69 


0.8 


1.2 


-21.85 


2.8 


-23.52 


0.7 


-20.69 


-22.06 


0.85 


1.1 


-21.98 


2.8 


-23.61 


0.7 


-21.11 


-22.49 


0.9 


1.1 


-22.11 


2.7 


-23.70 


0.7 


-21.79 


-23.18 


0.95 


1.1 


-22.24 


2.6 


-23.79 


0.7 


N/A 


N/A 


1. 


1.0 


-22.36 


2.6 


-23.88 


0.7 


N/A 


N/A 



a All five quasar surveys and lensing events are used. 



DM 



DM 



a 



DM 



203 km s" 1 , 
207 km s -1 , 
= 178 km s~ 



45/i 
28h 



pc, and 7 



al 

used 



30 , Maoz et al. P9| 



= 2.83. 
pc, and 7 = 1.89. 
= 45/i~ 1 pc, and 7 = 1.89. Based on surveys by Crampton et 
Surdej et al. 
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and Yee et al. ||3~2" |. Four lensing events are 
zi = 2.9157 in system Q1208+1011 is used. There is no solution for r c and <tdm of 
the system Q1208+1011 when fi G > 0.95. 
e zi = 1.66 in system H1413+117. 



Table 8: Expected lensing events and likelihood analysis. This table is for open universe 
models, i.e., = 0. N exp is our theoretical prediction for the number of lensing events, and 
ln£ is the logarithm of the maximum likelihood function. z\ = 1.438 in system H1413+117 
and z\ = 1.1349 in system Q1208+1011 are generally used, unless otherwise noted. 
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